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Abstract. In this paper we find explicit formula for some physically interest- 
ing solutions of multidimensional zero pressure gas dynamics system and its 
adhesion approximation. This includes the spherically symmetric solutions. 
First we write down solution of the adhesion model and study its asymptotic 
behaviour. The radial components of the velocity and density satisfy a simpler 
equation. This enables us to get explicit formula for different types of domains 
and study its asymptotic behaviour. A class of solutions for the inviscid system 
with conditions on the mass instead of conditions at origin also is analyzed. 



1. Introduction 

One of the analytical models proposed to describe the large-scale structure of 
the universe is the adhesion model 

ut + (u.V)u== |Au, pt-f V.(pu) =0. (1.1) 

Here u is the velocity and p, the density of the particles (see P] and references 
therein for the physical importance and analysis of solutions). This system and its 
inviscid counterpart, the multi-dimensional zero-pressure gas system 

Ut + {u.V)u = 0, pt+V.{pu) =0 (1.2) 

have been an active field of research since they were introduced, see jH 131 IH [171 IE] ■ 
However there is no clear analytical understanding of the solutions of these equa- 
tions. One question that remains is a well-posedness theory and large time be- 
haviour of solution. In this paper we construct special solutions which are physically 
interesting and study its asymptotic behaviour. 

The advantage of (jl.2|) over (jl.ip is that it is simpler to describe the solution 
in the region where classical smooth solutions exist. It is well-known that smooth 
solutions for (|1.2I) does not exist globally m {{x,t) : x £ > 0}, even when the 
initial data 

u{x,Q) = uo{x), p{x,0) = po{x) (1.3) 

are smooth with compact support. So the global solutions should be understood in 
a weak sense. Weak solutions are not unique and we need additional conditions to 
select the physical solution. One way of selecting the physical solution of (|1.2p is 
by taking the limit of solution of (|l.ip as e goes to 0. 

In the context of the large scale structure formation, the fastest growing mode in 
the linear theory has decaying vorticity. So it is natural to seek potential solutions 
of the model equations (11.11) and (|1.2p . Then the velocity can be represented in 
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terms of a velocity potential <j), see [18]. In this case, m*^, e > can be constructed 
using the standard Hopf-Cole transformation, see [8j [T^- Then the continuity 
equation for p'^ is a linear equation with smooth bounded coefficients which can be 
solved using the method of characteristics. We use this method to construct weak 
asymptotic solution of (II. ip and (II. 3p . 

For the inviscid system ()1.2|) . spherical solution of the form u = {x/r)q, p[x) = 
p{r), r — \x\, was constructed for n = 3 in ,12 . This is a special case of potential 
velocity. It was shown that the radial component of velocity and density satisfies 

qt + qqr =0, pt + \{r^pq)r = 0, r > 0, (1.4) 
where r = \x\. By the change of variable p = p, the equation (jl.4l) takes the form 
qt + qqr = 0, pt + {qp)r = 0, r > 0. (1.5) 

Classical theory of hyperbolic conservation laws of Lax pLS. is not applicable due to 
formation of 6 waves. However the system p.Sp has been well studied in [71 1^ 111). 
Using explicit solution of this system, we obtain radial solution of (|1.2p , with initial 
data at t = and with different behaviors at the origin x — 0. It is clear that 
in the region of smoothness, u and p constructed using (|1.5|) is a solution of the 
multidimensional system (jl.ip . As the solution is not smooth, it is not obvious that 
the constructed distribution is a weak solution of the multidimensional system (jl.2p . 
We show that along the surface of discontinuity, the Rankine Hugoniot condition 
derived by Albeverio and Shelkovich [1 is satisfied. 

In section 2, we construct weak asymptotic solution of (|1.2p . In section 3, we 
derive equation for the radial case of (jl.ll) and its linearization. In section 4, we 
solve boundary value problem for radial adhesion model and study its asymptotic 
behaviour for space dimensions n — 2 and n = 3. In the last section, we find explicit 
solution for the radial inviscid system, with conditions on the mass. 

2. Explicit weak asymptotic solution for multi-dimensional case 

In this section, we construct explicit weak asymptotic solution of the system 
(|1.2p with initial conditions (11.31) using the adhesion approximation (II. ip . with an 
additional condition that u = V xip- 

First we recall the definition of weak asymptotic solution from [T] who studied 
the system (|1.2p without the condition u — Wxtf)- A family of smooth functions 
(m'^,P'^)c>o is called a weak asymptotic solution of (|1.2p and (|1.3p if 

ut + {u\W)u^ = o'^(i?")(l), Pt + W.ip^u^) = o'^(i?")(l), 

u%x,Q)-uo{x)=o'^{R''){l), p%x,0) ~ po{x) ^ o'i,{l) 

where < o^(i?")(l), 77 >-j> as e -J> uniformly in i > 0, for every 77 G C^(i?"). 

We consider the adhesion approximation p.ip with initial conditions which are 
regularizations of (|1.3p . The velocity component u'^ can be constructed using Hopf- 
Cole transformation, and then the continuity equation for p'^ is solved using the 
method of characteristics. The resulting family (u*^, p'^)e>o, is shown to have neces- 
sary estimates to form weak asymptotic solution of l|1.2p and (|1.3p . 

Theorem 2.1. Assume uo{x) = where (po e W^'°°(i?") and po G L'^{R"). 

Let (f>Q — (pQ^rj'^ , V x4>h = x4>G*'n'^ md po = po*r]'^, where rf is the usual Friedrichs 
mollifier in the space variable x G i?". Further let 
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u'^{x, t) 



Vy0^(y)e ' L * \dy 



-1 i-y|- 



■+0o(y) 



(2.2) 



p^(x,t) = p;^(X^(a;,t,0))J^(a:,i,0), 

where X'^{x,t,s) is the solution of ^^^^ — u'(X, s) iwit/i X(s — i) = a; and 
J'^{x,t,0) is the Jacobian matrix of X'^(x,t^O) w.r.t. x. 

Then (u'^,p'^) is a weak asymptotic solution to (|1.2p and (jl.3p . 

Proof. Following [5j HJ [18] , we note that if 0'^ is a solution of 

0t + - Ia^/-, 0(x, 0) = 0^(x), (2.3) 

then u'^ = Va;0^ is a solution of with initial condition u'^{x,0) — Wx4io{x). 

Now using the Hopf-Cole transformation 61 = , it follows that (j>'^ is the solution 
of (1131) iff 9^ is the solution of 

9t^^Ad, dix,0)=e~^. (2.4) 

Solving (P^ . we get, 



Now 



^ ' (27rte)»/2 ^ ^ 



e^{x,t)x. = ^ \ / d.,-e-i^^+'^o(y'>^dy 
^ ' (27rie)"/2 _/^„ 



1 



a^^.(e-7^^1^).e-7<^o(y)dy (2.6) 



Vj 



(27rte)"/2 

-(^l^-'^°(^)^'''^-^"''^^^^'^^- 
In the last line we used integration by parts with respect to the variable y. Since 



= 6 ' , we get from p.Sp and (|2.6 



From this formula it is clear that 



k^lL°°(i?"x[0,oo)) < |Va;0oU°°(iJ")- (2-8) 

Further as differentiation under the integral sign is justified if \/x4'o is bounded 
measurable, it follows that u"^ is a C°° function in i?" x (0,oo). Indeed, we can 
write the formula ()2.7p for u*^ in the following form 



u'(x,t)= -"^ vv^w ^ 2.9 
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Differentiating p.9|) and using chain rule, it follows that 



Id^diu^l < (2.10) 

where Caj depends only on |Vx(/>o|l=°- Next we consider the continuity equation 
with coefficient it*^ and regularized initial condition, 

Pt + V,.KP) = 0, p%x,0) = p'oi^). (2.11) 

As u'^ is smooth, we use the method of characteristics to find p. Let X{x,t, s) 
be the solution of 

= u'{X, s) , X{s = t)^x. (2.12) 

as 

Since satisfies the estimates (|2.8p and (|2.10p . by the existence and uniqueness 
theory of ODE, there exists a unique solution X{s) to (|2.12p for all < s < <. To 
make the dependence of x and t explicit, let us denote it by X{x,t,s). This flow 
takes the point {x,t) to an initial point {X(x,t,0),0) and conversely. 

Let J{X'^(x,t, s)) be the Jacobian determinant of X'^{x,t,s) with respect to x. 
Then 

p'{x, t) = pl{X{x, t, 0))J{X'ix, t, 0)) (2.13) 

is the solution of (|2.1ip . The family {u'^,p'^)^yo given by (|2.7p and (|2.13p is a 
weak asymptotic solution. This follows easily as 

e / Au'^ri{x)dx = e / u*^ Ari{x)dx — 0(l)e 

uniformly in t for every rj € C^{R") as u"^ is bounded independent of e and 
{x,t) by the estimate (|2.8p . The solution (u^,p'^)e>o satisfies the initial conditions 
because we have 

u'{x, 0) - Wa^Mx) = u'{x, 0) - W^<j>'a{x) + V^(/)^(x) - V^Mx) 

= Vxipoi^) - Vxipoix), 

p'ix,0) - po{x) = p'{x,0) - p'q{x) +Po(a;) - po{x) = p^g{x) - po{x) 

where Vxf/^oix) — Vx^'oix), Po{x) — Poix) go to zero in distributions. 

□ 

There is a large class of interesting initial data which admit gradient type so- 
lutions as described in the theorem. One such case is when the initial data is of 
radial type. 

u{x, 0) = uq{x) = -qo{r), (2.14) 
r 

Clearly it can be written as a gradient, 

uo{x) ^ V x(l)o{x) ^ V qois)ds). 

Jo 

For density we give the initial condition 

pix,0) = Poix). (2.15) 
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Theorem 2.2. The solution of with initial conditions (|2.14p and (|2.15p with 

Jo Qois)ds < (X) and po{x)dx < oo /las the following asymptotic behaviour. The 
velocity component goes to as t tends to oo uniformly on compact subsets of R" 
and the mass p{x, t)dx is conserved. 

Proof. From ()2.5|) - (j2.6|) and = , we get the following formula for the velocity 
Making the change of variable y — )• (x — y)/V2t, we get, 



-1 Tl |2 , r|a=-\A2tT!/l / X , 1 

(. —\\V\ +fo ^ qo{s)ds\ 

u\x^ t) = ^-^^^ ^ (2.16) 

From ()2.16p . it follows that, as t goes to oo, u"^ goes to uniformly on compact 
subsets of i?" . From the formula (j2.13p we have 



p''{x,t)dx^ I p'Q{X''{x,t,0))J{X''{x,t,0)) 
poix)dx, 



since for every t, X'^{x, t, 0) : i?" — )■ i?" is a diffeomorphism and the determinant of 
Jacobian, J{X'^{x,t,0)) is positive. So the mass p{x,t)dx is conserved. □ 

3. Equations for radial case and their linearization 

To study the radial solutions of (jl.ip with prescribed initial and boundary con- 
ditions, we look for radial components of velocity and density of the form 

X 

u{x,t) = -q{r,t), p{x,t) ^ p{r,t), r = |a;|. (3.1) 
r 

We find the equations for q and p and show that the resulting system can be lin- 
earized using the Hopf-Cole transformation. More precisely, we have the following 

Theorem 3.1. The equation (jl.ip . for radial components can be linearized by the 
transformation 

q(r,t) = -e^, p(r,t)=r-("-i)p(r,i) (3.2) 
a 

to the following system 

€ ("- ^ 1) 

at = 7:["-rr-\ Cir], Pt + (pq)r = 0. (3.3) 

2 r 
Proof. A simple computation shows that 



^ ^ \{l/r)q + {xj/r^)qr~ ixj/r^)q, if, k = j 
^ \{xjXk/r'^)qr - {xjXklr^)q, if,k^j 



(3.4) 



and 



/y^ ^ f(^)9rr + (^-^)gr + (^-^)g, if,k^j 

\ 3)XkXk ] Xjxl ,Xj X,xl 2Xjxl X ,(3XJXI Xj . A f 1, / A ^' ' 
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From dsn]), ([331) and ([SH), we have for j = 1, 2, ...n 

n 

K)t + ^Uk{uj)^^ = {x,/r)qt + {xj/r)q[{l/r)q + (.T^/r'^)^^ - (x^/r^)^] 



k=l 



\ ^ Xk rXjX]^ XjX}^ 



—[it + qqr] 



and 



e ^ exjxj 3 3xj 3xj 3 , 

2 2 r r r-^ j.^' '> 



e . (n — 1) (n — 1) , 

= T,—\qrr + ^ -qr ~ —^q\- 

2 r r 
Using these in p.ip . we get 

Xj e Xj (n — 1) n — I 

—{qt + qqr\^7; — [qrr + qr :^q\, .?■ = 1,2,3,. ..n, (3.6) 

r Z r r r 

Also for the continuity equation, we have 

71 n 

pt + 5]Kp).. =Pt + Y.^{{l/r)q + {xl/r^)qr - {xl/r'')q)p + {xl/r^)qpr} 

k=l k=l 

= pt + {n/r)qp + qrp - {l/r)qp + qpr 
= pt + ((n - l)/r)qp + (qp)r 

(3.7) 

From (|3.6p and p.7p . it foUows that {u, p) of the form (j3.ip is a solution of (jl.ip iff 

q and p satisfy the equations 



gt + qqr = x 9rr H qr 2 — ^ ' 

2 r 



+ Z7iTT(-^""'V9). = 0. 



(3.8) 



Let p = r^" ^^p, then the above system, becomes 

e ('^ ^ 1) — 1) 1 

+ = 2 ilrr + , ^3 

Pt + {pq)r = 0, r > 0, i > 0. 

The equation for velocity component can be linearized using Hopf-Cole transfor- 
mation. For this first we note that if Q{r,t) is a solution of 

Qt + kiQrf)^^[Qrr+^-^^^Qr]. (3-10) 

2 2 r 

with initial and boundary conditions 

Q(r,0)= / qQ{s)ds, Qr{0,t)=qBit), 
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then q(r, t) — Qr{T, t) is a solution of the first equation of (|3.9p . with initial condition 
^(r, 0) = qo{r) and boundary condition q{Q,t) ~ qsit). Now let 

Q(r,t) = -eloga^ 

An easy calculation shows that Q satisfies (|3.10p iff a'^ satisfies 

e (n — 1) 

at = -[arr H Or , (3-11) 

z r 

□ 

We use these results in the next section to construct solutions of the radial system 
in domains with boundary. 

4. Radial solutions of adhesion model with boundary conditions 

Here we are interested in getting explicit solutions for in some spherically 
symmetric domains of the form fl — D x [0, oo) for 2 and 3 space dimensions. First 
we consider the domain ft = {{x,t) : \x\ — r < R,t > 0}. We take initial conditions, 

m(x,0) = i^goW, p{x,0) ^ po{r),\x\ < R, (4.f) 



and boundary conditions 



X 



lim —u{x,t) = qB, (4.2) 
\x\^B. \x\ 

where qs is constant. For p we need boundary condition only if qs is negative. 

lim p{x,t) = psit), ifqs < 0. (4.3) 

\x\-^R 

We also assume first order consistency condition of the initial and boundary data 
at {{x,t) ■.\x\=R,t = 0}. 

Let r = {{x,t) e fl : either = R or, t = 0}. Define pr on F by 

p,ir,t)^h^'^' ^^' = °' (4.4) 

As q is smooth, the characteristic of + {qp)r ~ passing through {r,t), j3{r,t,s) 
exists for all s < t, till it meets F, at time t — tQ{r, t) > 0. 

Let Jfc, for k — 0, 1, denote the Bessel functions of first kind and order k. With 
these notations, we have the following theorem. 

Theorem 4.1. Explicit solution for with initial and boundary conditions 

(liTTj) - is given by 

, X J,''drGiT,^,t)e-ifo'!o(s)'is^^ 

UiX t) — — 6 

l^^l J„'^G{r,tt)e-ifoiois)dsd^ (4 5) 

p{x, t) = r-("-i)pr(/3(r, i, t^{r, t))Mr, 
where G{r^ t) is described below. For n — 2, G has the form 

i-. CO n ^ A 2 

- ]^g((,.A)^i?^+.?jJo^(M.)''"^ir)^''(^)^ ^'-'^ 

where pn are positive solutions of the transcendental equation pJi{p) — kRJo{p) — 
and for n = 3, G takes the form 
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G{r,t,t) = — > ; — 77 r— r sm sin( )e 4.7 

i?r^^M?. + (teA)i?((9BA)i?-l) ^ ^ ^ ^ ^ ^ 

where fin cire positive solutions of the transcendental equation fiCOt{iJ,) + (qB /e)R — 
1 = 0. 

Proof. By Theorem 3.1, it follows that the velocity u is given by u{x, t) = — ejfy^, 
where a satisfies the linear problem 

e (n — 1) 
at = -\arr~\ aJ, r < R, t>0 

I r 

S,U(^)^s (4.8) 
a(r,0) ^ e - , r < R 

ear{R,t) + qBa{R,t) =0, t > 0. 

Then G given by ()4.6p and (|4.7p are just the Green's functions for the boundary 
value problem (|4.8p . for 7i = 2 and n = 3 respectively, see |16j. The formula (|4.5p 
for u follows. 

Now we find a formula for p. Notice that p{x, t) — where p satisfies 

qt + pqr = ~Prq- 

Method of characteristics says along the curve §7 = P, ^ = ^Prq- Let us consider 
the characteristic curve (/3(r, i, s), s) passing through (x, t). Integrating the equation 
for q, along this characteristic, we get q{x,t) = q{(3{r,t,to)) + J^^ ftoP'^P^^^^-"))'^'' ^ 
where io ~ toif, t){< t) is the time the curve touches the boundary T. □ 

Now we consider the equation (|l.ip in the domain ft = {{x,t) : i?i < |a;| < 
i?2, t > 0} with initial conditions 

u{x, 0) = ^(7o(r-), 0) = po(0, ^1 < l^^l < R2 (4.9) 

and boundary conditions 

lim ^u(x,t) = (7„z = l,2. (4.10) 

where qt are constants. For p we need boundary condition on (|a;| — Ri) only if 
is positive. 

lim p(x,t) =p,(t), i/, {-iy+\,>Q. (4.11) 

We also assume first order consistency condition of the initial and boundary data 
at {{x,t) : \x\ = Ri,R2a.ndt = 0} Let T = {{x,t) e Q, : either \x\ = R^ or, t = 0}. 
Define pr on F by 

Pr(r,t) = H';^' = (4.12) 
I^Pi(t), if r = 

As q is smooth, let /3(r, s) be the characteristic of pt + {qp)r — passing through 
(r, t) and let this characteristic meet the boundary point at time t = to{r,t). 
Let Jfe, k = 0,1, denote the Bessel functions of first kind and order k and Yk, 
k = 0,1, denote the Bessel functions of second kind and order k. With these 
notations, we have the following theorem. 
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Theorem 4.2. 



u{x, t) = — e 



r«2 



(4.15) 



(4.16) 



\A jl;^G{r,^,t)e-ifoioi-)'i^d^ (4.13) 
p{x, t) = r-("'i)pr(/3(r, to(r, t)), t^{r, 0)e-^4<..*) ^.(/^('^.M))^.^ 
where G{r,^,t) is described below. For n — 2,G has the form 

= YE#[(«2A)>/o(A„r2)-A„Ji(A„r2)]^ei^nM^r„(Oe-^ (4.14) 

n— 1 

S„ = (\l + (<Z2/e)')h(fci/e) Jo(A„i?i) + KJiiKRi)? 

- {\l + (fci/e)')[(g2/e)Jo(A„i?2) - A„Ji(A„i?2)]', 
Hn{r) = [-(gi/e)ro(A„i?i) + A„Yi(A„i?i)]Jo(A„r) 

- [~{qx/e)M\nRi) + A„ Ji(A„i?i)]yo(Anr) 
and the A„ are positive roots of the equation 

[-(gi/e)Jo(Ai?i) + AJi(Ai?i)][(g2/e)>o(Ai?2) - Ayi(Ai?2)] 
- [(g2/e)Jo(Ai?2) - AJi(Ai?2)]h(gi/e)>o(Ai?2) + AYi(Ai?i)] = 0. 
For n = 3, G has the form 
G{r,tt) 

^ (bi + fliA^Jvl/„(r)vl/„(e)e-4^^ 

r (i?2 - RiM + RlXDibl + i?2A2 ) + (6ii?2 + b2Ri){bib2 + R1R2XI) ' 

(4.17) 

w/iere 

*n(r) 61 sin[A„(r - i?i)] + i?iA„cos[A„(r - i?i)] 

61 = -(gi/e)i?i + l, 62 = (92/e)i?2 - 1, ^^'"^^^ 
and A„ are positive solutions of the transcendental equation 

(6162 - i?ii?2A^) sin[A(i?2 - Ri)] + X{Rib2 + i?26i) cos[A(i?2 - i?i)] = 0. (4.19) 

Proof. From the theorem (3.1), we get u — where a satisfies 

e (n — 1) 

at = 7:[arr + ^ -ar], i?i < r < i?2, t > 

z r 

SS^(=)^- (4.20) 

ear{Rt,t) + qia{R„t) ^ 0, i ^ 1,2, t > 0. 

Then G given by (14. 141) - (14. 161) and (I4.17p - (l4.19p are just the Green's functions for 
the boundary value problem ()4.20p for n = 2 and n — 3 respectively, see |16j . 
The formula (j4.13p for u follows. The formula for p is obtained by the method of 
characteristics the same way as in the previous theorem and details are omitted. □ 

Remarks on large time behaviour: 

It is easy to get the large time behaviour of the velocity from the formulas. First 
we consider the problem in the domain {{x,t) : \x\ < R,t > 0}. Then u given by 
(|4.5p - (|4.7p has the following asymptotic form : 



10 



PAL CHOUDHURY, JOSEPH AND SAHOO 



lim uix.t) — —e- — -— — 5 —. 

where for n ~ 2, 

fii being the first positive solution of the equation fiJi{ii) — fc_RJo(M) — 0. For 
n = 3, 

G(r,0 = -sui(— sm — , 
r R n 

fii being the first positive solution of the equation /icot(/i) + — — f = 0. 

Now consider the problem in the domain {{x,t) : Ri < \x\ < R2,t > 0}. Then u 
given by (|4.f 3p - ()4.f 9p has the following asymptotic form : 

rR2 o ^/ ^^ -i fS, 



lim u(x, t) — ^ , , „ , .. 

t^oo \x\ jl^^G{r,^)efoioi-)'i^dC 

where for n = 2, 

G(r,0-ei^i(0^i(0 
where Hi is given by (I4.15P with Ai being the first positive root of the equation 
(|4J6l) . For n = 3, 

G(r,0 = -*i(r)vl/i(0 
r 

where *i(r) = 6isin[Ai(r - Ri)] + i?iAi cos[A„(r - Ri)], bi = -{qi/e)Ri + 1, 
&2 = {<l2/^)R2 — 1, and Ai is the first positive solution of the equation (I4.f 8p . 

Asymptotic behaviour of p is more complex and depends on whether mass flows 
out or flows in. Consider the mass m{t) = Jjj p{x,t)dx. When p is spherically 
symmetric and D — {x : \x\ < R}, 



From the second equation of p.Sp . we have 

^ -q{R,t)p{R,t) 
and when p is spherically symmetric and D — Ri < \x\ < R2, 

- - {q{R2 , t)p{R2, t) ^q(Ru t)p{Ri ,t)). 
This shows that mass is conserved if velocity is zero on the space boundary. 
5. Spherically symmetric solution for the inviscid case 
The aim of this section is to find explicit global radial solution of 

Wt + (M.V)M = 0, pt + V.(pu) = 0, (5.1) 

with initial conditions 

-u{x,Q)^qo[r), p(a;,0) =poW =r-("-SoW, r = (5.2) 
r 

a condition on the mass 

p{x,t)dx = psit) (5-3) 
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and prescribed normal velocity at the origin, 

X 

lim —u{x,t) = qsit)- 



K-s-o r 



(5.4) 



For initial data uq{x), in the space of bounded variation, the vanishing viscosity 
limit of the velocity component u remains in the space of bounded variation, see 
[8]. The density p is generally a measure. In an important paper, local study of 
the Cauchy problem for (|5.ip and propagation of delta wave front was carried out 
by Albeverio and Shelkovich [Ij in the frame work of the weak asymptotic method. 
We summarize the definitions and results in [I] relevant to our work. 

A distribution (u, p) is called a generalized S- shock wave type solution if it is 
the distributional limit of a weak asymptotic solution {u'^,p'^), as e goes to 0. They 
started with an ansatz of the form 



u{x, t) = u+H{S{x, t)>0) + u^H{S{x, t) < 0), p = p{x, t) + e{t)5s{x,t)=o (5.5) 



where U-\-,u^ and p are smooth functions away from the surface S{x, t) = 0. They 
considered a smooth approximation {u'^,p'^) of (j5.5l) by a suitable regularization of 
the Heaviside function H and S measure. They showed that if this regularization is 
a weak asymptotic solution in the sense (j2.1l) . then {u, p) should satisfy the equation 
(|5.ip in the region of smoothness and along the surface of discontinuity S{x,t) = 0, 
the Rankine-Hugoniot condition, 



is satisfied. Here ug — ^ "~ ? being the velocities behind and ahead of the 

surface of discontinuity, respectively and Tt — {x : S{x,t) = 0}. 

Conversely, if there exists a distribution {u, p) and a .surface S{x,t) = satis- 
fying the above properties, then there exists a weak asymptotic solution [u'^ , p'^) of 
()5.1|) whose limit is {u,p). 

Further .starting from an initial discontinuous wave front satisfying an entropy 
condition, they constructed {u, p) and a surface S{x,t) = 0, satisfying (|5.1|) in 
the region of smoothness and the Rankine Hugoniot condition (j5.6p along surface 
of discontinuity S{x, t) = 0, using the work of Majda |14[ 115) on existence and 
stability of multidimensional shock fronts. This construction is only local. 

For the radial case we give a global construction of {u,p). We also do not 
require entropy condition on the initial data for u, as the equation is reduced to 
one dimensional Burgers equation and hence we use the well-known one dimensional 
theory. 

As a first step of construction of {u,p), we derive the equation for the radial 
case. An easy calculation as in section 3, shows that q and p ~ r^^-'^) p satisfy the 



St + us.^^S\r, = 0, — + Vr, {eus) - i\pu] - [p\us).V :,S\t, 



t •) 



(5.6) 



equations 



qt + qqr = 0, pt + [pq)r = 0, 



(5.7) 



with initial conditions 



q{x,Q) ^ q[){r), p{x,0) = po{r), r = \x\, 



(5.8) 



and boundary condition for q and integral condition on p: 




(5.9) 
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which is to be understood in a weak sense [2] ■ Weak formulation for the boundary 
condition for q is 

either 5(0+, i) = qsit) 
or g(0+,t) < and q^{Q+,t) < q^{t) 

and for p is 

if q{0+,t) > then a;„„i / p{r,t)dr ^ peit). (5.11) 



2,., - ^...2 (5-10) 



To describe the solution, we follow [5J [51 HH] ■ We introduce a class of paths in 
the quarter plane D = {{z,s) : z > 0, s > 0}. For each fixed {r,ro,t),r > 0, rg > 
0,t > 0, C{r,ro,t) denotes the following class of paths (3. Each path is connected 
from the point (ro,0) to {r,t) and is of the form z — (3{s), where /? is a piecewise 
linear function of maximum three lines. On C{r,rQ,t), we define a functional 

m {qBis^rds + \ f {^')ds. (5.12) 

We call f3o the straight line path connecting (rg , 0) and (r, t) which does not 
touch the space boundary a; = 0, namely {{0,t),t > 0}. Then let 

Air,ro,t) = J{M^ (5.13) 

Any P S C*(r, ro,t) — C{r,rQ,t) — /3o is made up of three pieces, namely lines 
connecting (rg, 0) to (0, ti) in the interior and (0, ti) to (0, on the boundary and 
(0, ^2) to (r, t) in the interior. For such curves, it can be easily seen from (|5.12p that 

/•*2 ('„„(' s'l + 'l^ 

m^J{r,ro,t,t,,h)^~J^ V ^'+2j7 + 2(r^- ^^'^^^ 

For curves P G C* (r, rg , t) made up of two straight lines with one piece lying on the 
boundary r = 0, we can write down a similar expression as in the earlier case. 
It was proved in [10] that there exists (3* G C* (r, rg , t) and corresponding 

ti{r,ro,t), t2{r,ro,t) so that 

B{r, ro,t)^ min{ J(/3) : /? € C* (r, rg, i)} 

= min{ J(r, rg, ^1,^2): < ti < t2 < t} (5.15) 
= J{r,ro,t,ti{r,ro,t), t2 (r, rg , t)) 
is Lipschitz continuous. Further 

m{r,ro,t) min{J(/3) : /3 G C(r,ro,t)} 
= min{yl(r, rg, i), B(r, rg, t)} 

and 

Q(r, t) = min{m(r, rg, t) + / gg(s)ds, < rg < 00} (5.17) 

Jo 

are Lipschitz continuous function in their variables. Further, minimum in ()5.17p 
is attained at some value of rg > 0, which depends on (r, i); we call it rg(r, i). If 
A{r,ro{r,t),t) < B{r,ro{r,t),t), 

g(r,0 ^^^' + f"^^''\o{s)ds (5.18) 
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and if A{r, ro{r, t), t) > B{r, ro(r, t, t), 

rroir.t) 

Q{x,t) = J{r,ro{r,t),t,ti{r,ro{r,t),t),t2{r,ro{r,t),t))+ / qo{s)ds. (5.19) 

Here and henceforth ro (r, t) is a minimizer in ()5.17|) and in the case of (|5.19|) , 
t2{r,t) = t2{r,ro{r,t),t) and ti{r,t) = ti{r,ro{r,t),t). With these notations, we 
have the foUowing resuh. 

Theorem 5.1. With ro{r,t), A{r,ro{r,t),t), B{r,ro{r,t),t),ti{r,t),t2{r,t) as de- 
fined above, define 

ri^I^, ^fA{r,ro{r,t),t)<B{r,roir,t),t), 

and 

p^^^^^ ^ i-Q,^,^Po{z)dz, tf A{r,ro{r,t),t) < B{r,ro{r,t),t), ^^^^^^ 
\uJn-iPB{t2{x,t)), if A{r,ro{r,t),t) > B{r,ro{r,t),t). 

and set 

Then the distribution (u{x,t), p(x,t)) given by (j5.20p - (j5.22p satisfies (jS.ip in the 
region of smoothness and (|5.6I) . along a discontinuity surface. Further it satisfies 
the initial conditions (|5.2p . mass conditions (15.31) and normal velocity at the origin 
([O)) in the weak sense ([5?5)) - (|5.1ip . 

Proof. Exphcit formula for the entropy weak solution q satisfying (|5.8p - (|5.10p is 
given in [6l [9] . This formula involves only a finite dimensional minimization, namely 
(|5.15p - (|5.17p in three variables {ti, t2, y). The formula takes the form 

rr-r^^ ifA(r,ro(r,t),i)<B(r,ro(r,i),0, 
\t=tT(^^ ifAir,roir,t),t)>B{r,ro{r,t),t). > 
To get the component p we consider the problem for 

P{r,t) = - p{s,t)ds (5.24) 

J r 

SO that 

p{r,t) ^dr{P{r,t)). (5.25) 
It is easy to see from (|5.7p . (|5.8p and (|5.11l) that P must satisfy 

Pt + qP^ = 0, 

with initial condition 

P{r,0)^- Pois)ds 

J r 

and boundary condition 

if q{0+,t) > then P{r,t)dr = ^—pB{t). 

An explicit formula for the solution of this problem is given in [5] and has the form 
(15.211) . The formula (|5.20p - (|5.22p is then obtained from the transformation ([XT]) 
and (jO^ - ljOSI) . 
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Now we show that {u, p) given by ()5.20p - ()5.22|) is a solution of ()5.ip in the 
region of smoothness. This fact follows easily as {q,p) satisfies the system ()5.7p 
and the weak form of initial and boundary conditions (|5.8|) - (I5.11|) (see 9 ). Also 
it was shown in 9 , that {q,p) is a weak solution of (|5.7p and that the surface of 
discontinuity of {u, p) and {q,p) are same. In order to show that {u, p) satisfies the 
Rankine-Hugoniot conditions, we use the corresponding conditions for the system 
(j5.7p . This follows from [l][T7] and we give the details here . 

Consider a solution (q{r,t),p{r,t)) of the system (j5.7p with a discontinuity on 
the surface S{r,t) = 0, with S{r,t) — r — s{t). In a neighbourhood of this surface, 
we assume {q,p) has the form 

q{r,t) = q{x,t), p{r,t) ^ p{r,t) + e{t)6r=s{t)- (5.26) 

where q,p are smooth except on the surface r — s{t), and e(t) differentiable function 
of t. The Rankine-Hugoniot condition for {q,p) takes the form 

ds{t) q+ + q^ de ds 

We show that the distributions ([Og)) - ([0^ satisfy As u = f g, A"--'^^ p = p, 

the coefficient of the measure S in the density p, e and that in p are related by 

r("-i)e = e. (5.28) 

Since S(x,t) = S{r,t) ~ r — sit), with r — \x\, we have 

dt r 



The first equation of (15. 6p becomes 

e 1 Y7 CI ds q+ + q^ X X ds q+ + q^ 

+ us.V.Slr, = + = + - 0, (5.29) 

where we used the first equation of (|5.27p .To verify the second equation in (j5.6l) . 
we compute each terms in radial components. 

Se _ de de _ de St de _ de 

St^dt^^dr^dt Srdr~ df ^^-^"^ 
Also we know from [T], 

VT,.{eus) = -2KGe, (5.31) 
where K is the mean curvature of the surface of discontinuity, K = — ^V.i^, 



V 



and G = — -i^. An easy computation gives 



22^dx,~ 2r^ 2^r^~ 2r ' ^^"^ ' 

Z— 1 2—1 

Using ([02)) in ([QT]) . we get 



n — 1 , St , n — 1 dr 
Vr..(e..) = ^(-^)e.^-e. (5.33) 



Finally, 



(5.34) 
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Since 5'^ = 1, we get from (lOOl) . (|533l) and (fS^i]) . 

— + Vr,(eu5) - (H - [p]u5).V,^|r, = ^ + ^e- - [p] — ) 

- [pg] - [p] 



where in the last equality we used the second equation of (I5.27p . This completes 
the proof of the theorem. □ 
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